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Abstract 

We make a comprehensive analysis of an extended supersymmetric model(ESSM) obtained by 
adding a pair of vector-like families to the minimal supersymmetric standard model and having 
specific forms of 5 x 5 fermion mass matrices. The singlet Higgs couplings which link the ordinary to 
vector-like generations do not have the renormalization effects of the gauge interactions and hence 
the "quasi-infrared fixed point" near the scale of the top quark mass. The two-loop Yukawa effects 
on gauge couplings lead to an unified coupling ax around 0.2 with an unification scale Mx of 10^^'^ 
GeV. Large Yukawa effects in the high energy region arrest the growth of the QCD coupling near 
Mx making the evolution flat. The renormalization effects of the vector-like generations on soft 
mass parameters has important effects on the charge and color breaking(CCB)minima. We will 
show that there exists parameter space where there is no charge and color breaking. We will also 
demonstrate that there exists minima of the Higgs potential which satisfies the mass of the Z boson 
but avoid CCB. Upper limits on the mass of the lightest Higgs boson from the one-loop effective 
scalar potential is obtained for sets of universal soft supersymmetry breaking mass parameters. 
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1 Introduction 



The commonly accepted notion of a unified coupling ax = 0.04 at 10^^ GeV has been 
questioned recently in view of a conjecture that the four dimensional unified string coupling 
is likely to have a semi-perturbative value (0.2-0.3) at the string unification scale (see later) 
so that it may be large enough to stabilize the dilaton but not as large as to ruin the 
supersymmetric gauge coupling unification [Q]. 

It may be noted that the MSSM unification scale Mx = 2 x 10^^ GeV appears to be 
a factor of 20 smaller than the one-loop level string unification scale of Mst ~ Qst 5.2 x 
10^'^ GeV ~ 3.6 X 10^^ GeV[^ In this context the extensions of MSSM spectrum addressing 
the issues of ax and Mx may be interesting. The specific extension which we will consider 
is called as the Extended Supersymmetric Standard Model (ESSM). This itinerary extends 
the MSSM particle spectrum by adding two light vector-like families which are Ql,r = 
{U, D, N, E)l^r and Q'^ ji = {U' , D' , N' , E')l^^ and two Higgs singlet scalars Hs and Hx- 
These new particles and as well as their superpartners are 1 to a few TeVs[|l[. ESSM has no 
unified gauge group as such and so the corresponding leptoquark and di-quark gauge bosons 
leading to proton decay does not exist. The extra sets of matter fields denoted by previous 
authors Ql\Q'r and Qr\Q'l can be thought of as 16 and 16 representations of a cosmetic 
unified gauge group of SO(IO). 

It has been noted sometime ago that precision measurements of the oblique electroweak 
parameters and the number of light neutrino species does not favor extra chiral fami- 
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lies. They do not however constrain vector-hke famihes due to decouphng effects|Q, ^. Such 
vector-hke famihes i.e., 16 and 16, which apparently are predicted genericahy in string the- 
ories could well exist in the TeV region. It may be worth understanding that the pair of 
vector-like families could give us a clue to the family mass hierarchy within the three chiral 
generations via a see-saw mechanism residing inside the 5x5 fermion mass matricesQ. 

Such a pair of vector-like families and thus ESSM provide some unique advantages over 
MSSM. They are as follows. 

• One-loop evolution of three gauge couplings maintains the approximate meeting when 
complete extra families such as vector-like families are included. 

• Although one-loop approximation leads to the unification of the three a^'s at the same 
scale Mx the unified coupling ax increases. 

• Numerically speaking ax can be raised to 0.2 — 0.3 in one loop. Two-loop gauge couplings 
accelerate the growth whereas the two loop Yukawa effects on gauge couplings does the 
contrary. This is vital for having a good fit to as{mz) within unification. 

Now we briefly describe the work of the previous authors[|l|. Instead of using the 6 
function approximation to study the particle thresholds they used exact threshold functions 
[]^, ^. The evolution of the Yukawa couplings was studied in one-loop. The off-diagonal 
Yukawa couplings between the third family and the vector-like families and the self couplings 
of the vector-like families were taken to be large varying between 1.0-2.0. This was done for 
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the Yukawa couplings of the up, down, charged lepton and neutrinos so that they approach 
their "quasi infrared fixed point" values near the scale of the top quark mass. The combined 
effects of the two-loop gauge contributions with one-loop Yukawa contributions and threshold 
effects exhibited the following features. 

• Three couplings still met. This was for a wide range of variation of the particle spectrum. 
A higher ax ~ 0.2 — 0.3 as well as higher Mx ~ 10^'' GeV compared to the MSSM values 
ax ~ 0.04 and Mx ~ 2 x 10^^ GeV was found. ESSM unification scale Mx is now closer 
to the string scale. The remaining gap of a factor of 5 or so between Mx and M^t may 
not be alarming. It could partially be due to the increased influence of the two-loop string 
threshold effects. This is because a semi-perturbative range of values of ax is at work. Such 
a value can bring a correction to the one-loop formula of Mgt- 

• The couplings meet for lower values of Ciz{mz)\jfs^^^^ ~ 0.120 — 0.124 in the ESSM case. 
This is in good agreement with the experimentally observed value of az{f^z)\°MS ~ 0.117 -|- 
±0.005. By contrast we typically need higher values of a^{mz) such as a2,{rnz)\^^^^^ > 0.125 
if rrig < 1 TeV and M1/2 < 500 GeV in the MSSM case If. 

These promising features serve the motivation to study ESSM further. The purpose of 
this paper is to explore the unification of ESSM gauge couplings further by including two- 
loop Yukawa evolution(see figure l.b). Three- loop evolution of the three gauge couplings 
including vector-like families with contributions only from gauge interactions has already 
been studied to some extent 0. Our reason to study the two-loop evolution of the Yukawa 
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couplings of ESSM is that there are 20 Yukawa couphng entries. These couphngs can be 
near the perturbative maximum Q at Mx- We also expect that the two- loop contributions 
of the semi-perturbative gauge couplings to the evolution of the Yukawa couplings would be 
significant | [ill . 

It seems to us that it is important to ask the question whether the inclusion of the 
two-loop contributions of the gauge couplings to the evolution of the Yukawa couplings and 
vice- versa would preserve, improve or adversely alter the features of "precision" coupling 
unification studied previously |jl]]. Thus we think that two-loop gauge evolution coupled to 
two-loop Yukawa evolution is an interesting approximation to explore in the context of semi- 
perturbative unification. 

2 Evolution of couplings at two-loops 
2.1 Definitions 

We perform a two-loop analysis of the renormalization group equations (RGE) of the gauge 
couplings including the effects of the third generation Yukawa couplings at two-loops and 
carefully taking into account the threshold effects due to the spread of the superpartner 
masses and the masses of the particles of the vector-like families. Furthermore we use mass 

dependent running of the gauge couplings as given in(0) or equivalently stated as the effective 

^ This may be viewed in comparison to the non-perturbative values of the unified couphngs px| ] . 

^ We would estimate the ratio of two- loop versus one- loop contributions to the (3 functions for the Yukawa 
couplings to be ~ [{g^ orh?) / A-k\{Ati)~^ C, where the coefficients C are typically large. Taking a representative 
value of C ^ 10 — 30 it suggests that the ratio of relative contributions can be about 25 to 75 % near Mx 
for g747r - h'^/An - 0.25. 
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coupling formalism which guarantees a smooth cross over of the beta function coefficients at 
the threshold of each individual particle. The two-loop RGEs for the Yukawa couplings can 
be found in the Appendix and those for the gauge couplings can be expressed as 

7-1 7 T (fc) 

We have redefined t = In/x and Yk = hi/ An. The coefficients bi, bij and are given in 
the Appendix. In a mass dependent subtraction procedure these coefficients are replaced 
by threshold functions^. The threshold functions depend on the scale through the ratios 
ntk/ where is the mass of the particle running inside the loop. Integrating (|l]) we obtain 
the analytical expression of the effective couplings. At one-loop approximation they are 
given by 

a-\^^) = a-\mz) + E ^ (i^fc(mfc^) - Fu{mk/mz)) . (2) 

The threshold function Fk{mk/ ^i) has the value \a[mk/^j) in the limit rrik/fi 0. Therefore 
when all the masses are well below the scale n we recover the familiar one-loop expression 

a-\f,) = a-\mz)-^\n(-^). (3) 

ZTT \mz/ 

On the other hand Fkijnk/ jj) vanishes in the limit mk/n — > oo showing the familiar decou- 
pling of heavy masses from the running of the couplings. A smooth threshold crossing with 
the variation of /i is obtained thereby. In our approximation the two-loop parts of the beta 
functions bij and changes at each threshold in a step-function approximation though. The 
details of our procedure to determine the masses of the superpartners with the constraint of 
correct radiative electroweak breaking in ESSM will be given in the following sections. 
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At the unification scale Mx we have the Yukawa couphng matrix of the third and the 
heavy vector-hke generations 16 and 16 in the simple form 



I63 16 16 



I63 



He - 16 



16 



(4) 



XfHf yfHs 
y'^Hs ZfH^ 

\x',Hf z',Hx J 
The couplings Xf, x'^ {yf, y'J denote the interactions between the chiral and the vector- 
like generations via the doublet scalars Hi, H2 and singlet scalar Hg whereas Zf, z'^ give 
the interactions within the vector-hke generations through the singlet B.\. Note that this 
is only a formal representation of the Yukawa matrix. In practice they should be viewed 
in terms of the component Yukawa matrices given in the Appendix. When we run down 
the couphngs from the unification scale to the electroweak scale the vector-hke generations 
become massive and consequently the Yukawa matrices get rotated projecting the effective 
third generation Yukawa couplings ht, hb and hr at niz- Hence the RGE for the Yukawa 
couplings are integrated in a top-down approach in the ranges Mx — > Mq — > M^ — > mz, 
with the boundary conditions at the vector-like scale 



hr{ML 





+ 












+ 








Zq 




+ 





Mr 



Mq 



(5) 



Mi 



Here Mq denote the quark and Ml denote the leptonic members of the extra generations 
consequently. Furthermore we call vs =< Hs > and v\ =< Hx >. Assuming all the 
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Yukawa couplings to be large at the unification scale {hi{Mx) = -\/47r) we evolve them to the 
scale mz- For the sake of a compact notation let us give it the name complete fixed point 
scenario(CFPS). We fix Vs/v\ ~ 0.5 for the time being. This prevents a low value of ht{mz) 
in terms of the ESSM couplings described in (j^) and consequently prevents a low prediction 
of the top quark mass. Later we will prove that there exists charge and color conserving 
minima with our choice of vs/v\ ~ 0.5. Furthermore we see that the singlet Higgs couplings 
do not feel the QCD interaction and thus cannot approach a "quasi-infrared fixed point" 
near the scale of the top quark mass and thus we may suspect that the radiative electroweak 
breaking may occur at the wrong place if their initial value at Mx is taken at \/47r. 

Imposing the unification condition ai(Mx) = a2{Mx) = a^{Mx) = ax the gauge cou- 
plings are evolved down to niz scale. We use the central values of sin^ 9^ = 0.2319 and 
Ci~^{mz) = 127.9 to fit Mx and ax- Now Mx and ax are fitted (which depend on the mass 
spectrum of the vector-like generations and superpartners as well as Yukawa couplings) so 
we predict the QCD coupling asi^niz) by running back using two-loop RGE. 

To compare the values of ai{mz) with the experimental ones we have consistently trans- 
lated the effective couplings to the values in the MS scheme. We have found that in the 
present case U{1) and SU{2) effective couplings are only 1% larger than the ones in the MS 
scheme. However a^{mz) can be 5-8% larger depending on the supersymmetry spectrum 
under consideration. 
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Figure 1: (l.a) Running of the gauge couplings when the Yukawa couphngs are taken at 
one loop (dash line) and at two-loop (solid line), (l.b) Running of the Yukawa couplings at 
one-loop (dash line) and two- loops (solid line). M2 = 90 GeV, rhg = rfih = 1 TeV, Mq = 3 
TeV. 

2.2 Comparison of one and two- loop results 



Let us choose a feasible superparticle mass spectrum which can be used to evaluate the 
threshold functions. Later we will calculate the evolution of the soft mass spectrum more ac- 
curately using universality and use iterative numerical subroutines to calculate the threshold 
functions. We choose Mq — 3 TeV, wino mass M2 — 90 GeV, squark mass m^^ and higgsino 
mass rhh as 1 TeV. The more massive Higgs scalars are also taken at 1 TeV while keeping 
the lightest Higgs mass at 95 GeV. We have neglected the mixing between the charginos 
and higgsinos as well as between left and right handed squarks and sleptons which occur in 
the threshold functions. We have assumed the mass pattern of the vector-like generations as 
Mq ~ niQ and Ml ^ ttil- This reflects the assumption that supersymmetry breaking scale 
is near 1 TeV whereas the vector-hke scale is few TeVs. Thus we are led to the approximate 
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relation Mq ~ SM^. The factor of 3 within the vector-hke generations represent the similar 
QCD renormalization effects parallel to the three chiral generations. In figure ^a we have 
compared the evolution of the effective gauge couplings coupled with the Yukawa couplings. 
Yukawa couplings are calculated at one-loop(dash line) and two-loops(solid hne). 

We first notice that in the high energy region the two-loop Yukawa couplings are ap- 
preciably different from the one-loop ones (figure [^.b). The curves for the gauge couplings 
are flattened near the unification scale. Note that effect is pronounced in the case of the 
QCD coupling. Due to the presence of a large number of Yukawa couplings their contribu- 
tions to gauge coupling beta function coefficients becomes comparable the gauge coupling 
contributions. They may even be mutually canceling causing a flatness of the curve. This 
cancellation is more forceful when the Yukawa couplings are integrated at two-loops. At low 
energy the Yukawa couplings approach their "quasi-infrared fixed points". The distinction 
between the one-loop and the two-loop results for the Yukawa couplings reduce in the low 
energy region below the vector-like threshold. From the point of view of gauge coupling 
unification the overall effect would be a small increase of the values of and 0:3 ^(m^) and 
Mx. 

Stage is set for the question, "how much these predictions vary when the scale Mq 
is varied ?" Intuitively it is clear that we should recover the MSSM results in the limit 
Mq — ^ Mx- In figure |^ the predictions of ax and a^{mz) are plotted when varying Mq with 
Yukawa couplings at one- loop (curves labelled (a)) and two- loop Yukawa couplings (curves 
labeled (b)). They also describe the results including the threshold effects due to a higher 
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Figure 2: The predicted values of a^i^niz) and ax with varying Mq. The sohd and dash 
hues are for different M2 as given in the legend. Curves labelled (a) are obtained with 
the Yukawa couplings integrated at one-loop whereas those labelled (b) shows the Yukawa 
couplings at two-loops. 

gaugino mass of M2 = 200 GeV. Once m^^ and M2 are fixed so are Mq and M1/2 at the 
unification scale (figure 

The plotted values of as^mz) are their MS scheme value. When Mq is near the TeV range 
the three gauge couplings are non-asymptotically free for most of the range of integration 
and later they merge at a larger value of ax compared to the MSSM prediction. When Mq 
increases asymptotic non-free behavior reduces and the MSSM results are recovered. 

2.3 Fixed point scenario, vector-like scale and fermion masses 

Now we have to choose an input value of the ratio vs/vx to translate Yukawa couplings to 
fermion masses using (^. We continue with our choice of the ratio vs/v\ ~ 0.5. Next we 
have to fix all Yukawa couplings of isospin up sector {H2 coupling) as well as isospin down 
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Figure 3: The figure at left gives tlie predicted values of Mx', labels (a) and (b) imply Yukawa 
couplings integrated at one-loop or two-loops. The figure in the right gives the extrapolated 
values of Mq and M1/2 (at the unification scale) requiring that the lightest charging mass 
Ml/2 is 90 GeV. First and second generation squark masses are consistently taken at 1 TeV. 



sector (i^i couphng) at Mx- This is done by choosing both sets at v47r. After evolving the 
Yukawa couplings the related fermion mass scales we input uIt- = cos/3 ^ = 1.777 GeV 
for the isospin down sector and calculate the famihar tan /3 which is of MSSM. At this 

stage we can check the predictions of the top quark mass and the bottom quark mass. We 
have found mf'^ = 165 - 180 GeV, mb(mb) = 4.1 - 4.4 GeV. The range corresponds to a 
wide variation in Mq is in the domain 3 TeV-10^^ GeV. Each value of Mq gives a value of 
tan /3 via rrir- The range of tan 35 to 60 is simply reflecting that all the Yukawa couplings 



at Mx are v47r. 



An interesting experimental situation from the point of view of semi-perturbative uni- 
fication would be to have the vector-like scale of a few TeVs. The mass of the vector-like 



particles can be related to an "effective" parameter. We will define it in a moment. A sim- 



12 



pie numerical estimation will show us that if we demand that the vector-like squark masses 
are at the TeV region the scenario with all the ESSM Yukawa couplings are \/47r at Mx is 
excluded. It will be apparent that the reason is related to the link between ^ and v\. 

The superpotential for the Higgs scalars is 

W = hH,H2H^ + jHl + ^Hl (6) 

This specific superpotential is an inspired guess as suggested in|]l|. The scalar Hs plays little 
role in the EWSB as it does not couple to Hi and H2. However the singlet Hx plays the role of 



the field N in the Next-to-Minimal Supersymmetric Standard Model (NMSSM) [0, ^ |14 
and EWSB is driven when gets a VEV (see section 3.2 and figure 4 later). An effective /i 
parameter can be defined by the relation 

fi = kivx . (7) 

The coupling ki is computed near mz scale. We expect that ~ 1 TeV. The vector-like 
fermion masses are obtained as the eigenvalues of the matrix given in (^). In the first 
approximation they are 

and similarly for the heavy leptons, where the ratio z'^{Mq) / z[{Ml) ~ 3 is found due to 
renormalization effects. It is very easy to see that we have the "quasi-infrared fixed points" 
after integrating the RGE ^ 

4(Mq)~ 0.73, ;2„(Mq) ~ 0.56, Zd(Mq) ~ 0.53, and fci(MQ) ~ 0.01. (9) 

^ ki does not feel QCD interaction however it still feels the SU(2) interaction via Hi and H2 which cannot 
pull it up too much when the scale approaches mz scale from above. 
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We see vx = ~ 100 TeV and thus the vector-hke particles are at a scale ~ 50 TeV. This 
is far too large to be interesting. Therefore if we consider smaller values of the vector-like 
scale the initial values of the couplings Zf and z'^ at the unification scale have to be smaller 
so as to fit the input value of Mq of few TeVs. This is our first reason to reduce the Yukawa 
couplings from the CFPS scenario. As a matter of fact smaller Yukawa couplings are forced 
by potential minimization conditions for any value of Mq as described later. 

In the next section we will study the soft mass evolution. But beforehand we may take 
note that the effect of large Yukawa couplings is to diminish the superparticle masses. Large 
Yukawa couplings will dictate charge and color breaking minima driven by cubic A couplings 



even when the weak constraint [O, IIH] 



\Aijk\^ < 3{m^i + m] + ml), (10) 

is satisfied. We must stress that especially when all the Yukawa couplings are at the "quasi- 
infrared fixed point" region near the scale of the top quark mass a detailed evaluation of the 



minimization conditions are necessary ||16|| . We will demonstrate in this case there is a charge 
and color breaking minimum deeper than the electroweak symmetry breaking one in ESSM 
for any value of Mq including the MSSM limit. Consequently we will conclude that one or 
more Yukawa coupling must be reduced at the unification point to reach a global charge and 
color conserving minimum. 

3 Soft mass spectrum of ESSM 
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3.1 Soft mass evolution 



Assuming the universality of soft supersymmetry breaking terms at the unification scale the 
supersymmetry spectrum at the weak scale is given in terms of following parameters at the 
unification scale apart from the dimensionless gauge and Yukawa couplings 

gaugino mass = M1/2 ; scalar mass = Mq ; cubic coupling = Aq. (11) 

The gaugino masses Mi are calculated integrating the two-loop RGE 

dt -^Sn^^U.- ^ ' 

The squark and slepton masses of the first and second generations are calculated neglecting 
small D-term contributions from the one-loop expression of 

drhl 1 



-—Y.AC,iRa)a,Ml (13) 



dln/i 271 

We have C2{Ra) = {N'^ — 1)/2N for the fundamental representation of SU{N). When 
calculating the squark and slepton masses of the third generation we have properly taken 
into account the effects of the related Yukawa couplings. Here we take as the initial 
value at the unification scale for those Yukawa couplings of ESSM which are related to ht{mt) 
below the vector-like threshold. This fixes the top quark mass. 



The one-loop RGE as in (|I^) can be easily integrated numerically. It is convenient to 
represent them after integration by the expression 

ml = Mi + CaMl/^ . (14) 
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The constants Ca depend on the values of the gauge couphngs and their /^-functions and 



thus imphcitly(via threshold effects) on the vector-like scale Mq. We again rewrite (|1^) in 
terms of the wino mass M2 instead of the universal mass parameter M1/2 in the following 
expression 

Thl = M^ + CaM2. (15) 

The values of Ca are much larger than those of MSSM. Sample values of Ca are given in table 
§ The MSSM limit is Mq ~ 10^^'^ GeV is also quoted in table for instant comparison. 
ESSM spectra gives an enhancement of the squark and slepton masses for a given M1/2 
compared to MSSM. This was noticed previously!^]. 

Table 1 gives such coefficients for the masses of the left handed squarks , left and right 
handed sleptons. We also give the coefficient for the ratios of the gaugino masses. Here we 
have taken the initial values Yukawa couplings at the largest allowed value (see illustration 
after (|^)). Of course later on we will fix the Yukawa couplings via potential minimization. 
For low values of Mq the squark and slepton masses of the first and second generations 
get large gaugino contribution due to the absence of the reverse effects from the Yukawa 
couplings. Hence even if we assume that the wino is as light as the experimental lower limit 
pushing Ml/2 to the least, the squark masses would still move to the TeV range . The left- 
handed slepton mass is larger than the gluino mass M3 and the right-handed slepton mass 
is roughly of the same magnitude as M3. As Mq increases we reach the MSSM limit where 

light gaugino masses together with large squark masses imply a large value of Mq as can be 

Those for the right handed squarks, up and down, are roughly the same as cq. In the numerical 
calculations of the effective couplings we have assumed degenerate masses of left and right handed squarks 



where both are taken to be ihq. 
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Mq 




Cl 


Cr 


M1/2/M2 


M3/M2 


M2/M1 


3TeV 


110.2 


21.2 


8.8 


8.8 


3.8 


1.6 


300 TeV 


27.3 


4.7 


1.8 


4.0 


3.9 


1.8 


10^0 GeV 


10.7 


1.2 


0.4 


1.9 


4.0 


1.9 


10^6-5 GeV 


9.3 


0.8 


0.2 


1.2 


4.0 


1.9 


MSSM 


9.3 


0.8 


0.2 


1.2 


4.0 


1.9 



Table 1: Coefficients for the left-handed squarks, left and right handed sleptons masses 
of the ffist and second generation. The MSSM limit is at Mq ~ 10^^'^ GeV. 

seen in figure |^. It is interesting to observe that the ratio ~ 4 is practically independent 

of Mq. 



3.2 Radiative electroweak breaking 



The Higgs sector of the ESSM resembles to that of the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM). The VEV of H\ generates effective fi and B parameters. In 
other words at low energy we evaluate the expressions 

fi = kivx, (16) 
B = hvx + A,. (17) 

Ai and A2 are the cubic dimension-full couplings corresponding to the dimensionless Yukawa 
couplings ki and k2 of the superpotential W. 

In MSSM there are two minimization conditions. They are the equations obtained by 
minimizing the scalar potential with respect to the real components of Hi and H2. Using 



these conditions one can fit two parameters and B with two inputs tan/3 and m^|jT^. In 
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Figure 4: Fitted values of the effective parameters fi = ki{mz)v\ and B = k2{mz)v\ + Ai at 
Low energy. For each parameter the lower curve correspond to Mq = and the upper curve 
to Mo = 800 GeV. We have taken kl{Mx) = 47r. 

ESSM there are three minimization conditions. The conditions for the minimization with 
respect to the real components of Hi and H2 along with the extra condition 



sin 2/3 = 2vx 



mjj^ + k2A2Vx + kfv^ + 2k^vl 
kiAiv"^ + 2kik2v'^v\ 



These three conditions can be used to fit three parameters (/c2, ^a) with three inputs 



(tan/5, fci(Mx), "m^). We have taken ki{Mx) = v47r. Once Aq is fixed the effective jj, and 
B parameters were computed using ([T6|) and (0) for various values of Mq keeping M2 = 90 



GeV. The values are displayed in figure ^. The value of /x is roughly constant with the 
variation of tan j3 at the TeV scale. The B parameter decreases with tan 13. 
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3.3 Superparticle masses of ESSM 



Let us remind ourselves that while calculating the unification coupling and unification scale 
(in the figures ^ and ^ ) we did not take into account the constraint from electroweak 
symmetry breaking (EWSB). We showed how the predictions depend on the vector- like scale. 
All the Yukawa couplings were taken as (v^ivr) at the unification scale. We see from(P) that 
in such a scenario Mq becomes too large to be experimentally interesting. In this section 
we will go beyond this approximation and move away from the fixed point scenario up to 
the degree to which correct electroweak breaking is obtained. This is our initial reason to 
reduce the values of the couplings accordingly. Later we will see that such a reduction has 
stronger basis as it is welcome from the point of view of charge and color conservation. 

In addition to this the complete fixed point scenario always leads to a large value of tan (3. 
If we want to consider ranges of values of tan (3 the initial values for the Yukawa couplings 
related to the bottom-tau sector (x/, x'^) should be tuned. This causes difficulty with our 
numerical subroutines. We will rather bypass this by taking a bottom-up approach taking 
tan /3 as a free input parameter at low energy. The initial Yukawa couplings are taken such 
that they reproduce the experimental values of mb{mb) = 4.4 GeV and = L777 GeV. 

Now we will describe our procedure of pinning down the Yukawa couplings for the EWSB 
mechanism. In table ^ we present sets of values of tan j3 and the corresponding fitted values 
of the Yukawa couplings. We also give the solutions for as, Mx and ax- While fitting we 
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tan /5 = 5 


tan/3 = 30 


tan/3 = 50 


v\/vs 


0.325 


0.317 


0.255 


Xd{Mx), Xr{Mx) 


0.013, 0.029 


0.096, 0.212 


0.449, 0.908 


ZgiMx), Zg{MQ) 


0.133, 0.351 


0.131, 0.343 


0.086, 0.253 


Zu{Mx), Zu{Mq) 


0.355, 0.621 


0.338, 0.604 


0.207, 0.469 


Zd{Mx), zaiMg) 


0.218, 0.621 


0.214, 0.604 


0.155, 0.469 


k2{Mx) 


0.932 


0.940 


0.297 


ki{mz), k2{mz) 


0.242, 0.158 


0.220, 0.172 


0.149, 0.113 


S3 


0.124 


0.123 


0.120 


Mx 


16.94 


16.93 


16.88 


ax 


0.252 


0.247 


0.226 



Table 2: Values of Yukawa couplings and the ratio v\/vs satisfying the fermion masses. We 
have taken Mq = 3 TeV, rrit = 175 GeV, rrib = 4.4 GeV, rrir = 1.7 GeV. The predictions for 
as, Mx and ax are also included. EWSB is guaranteed for these sets of Yukawa couplings. 
These values will be used to get the superparticle spectrum given in table || 

have chosen the following boundary conditions at the unification scale 

Xf = x'f , yf = zf, 

Zg = z[ , Zu = Zu, Zd = Ze. (19) 

The remaining Yukawa couplings are taken at y/An. Fitted values of a2,{mz) and ax de- 
crease with tan/?. The actual values can be found in table |^. In the large tan/3 region the 
Yukawa couplings related to the bottom-tau sector are larger and we get the limit of all 
Yukawa couplings in the fixed point region (section 2). Thus we have found sets of Yukawa 
couplings(table giving the correct EWSB and mz and shows that the mechanism works 
for ESSM as well. 



Now we describe the supersymmetry spectrum of ESSM including EWSB. The results 
are reported in table |^. We have calculated the masses of the Higgs scalars including the 
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mixing between the doublets and the singlet Hx. We have calculated the masses of squarks 
and sleptons, gauginos (chargino neutralino gluino) while guaranteeing a correct EWSB. In 
the first row we give the mass parameters required to obtain the spectrum. The constraint of 
correct EWSB leads to an effective /i parameter in the TeV range. We note that in the case 
of ~ 1 TeV and low values of M2 and Mi the gaugino-higgsino mixing is almost negligible. 
The heavier chargino and neutralinos are almost higgsinos. The heaviest neutralino is almost 
a singlino. The spectrum is not very different from that with an assumption of neglecting 
the mixing and requiring a higgsino mass of few TeVs. 

The same feature of small mixing is present in the Higgs sector as well between the 
doublets and the singlet. There are three CP-even states 5*1, 5*2, 5*3, two CP-odd states ai, 02 
and two charged states H^. Due to the large value of /i the Higgs masses coming mainly from 
the doublets, Si, S2,ai, can be thought of as h, H, ai, Higgs scalars of the standard 



MSSM Higgs spectrum|[T8|l. We have found that the predominantly doublet Higgs scalars 
satisfy the mass relation 

"^51 < rns2 ~ ruai ~ mH± ■ (20) 

Notice that the lightest (tree-level) Higgs mass quoted in table ^ is too low to be com- 
patible with the current experimental lower bound of 95 GeV. Higgs masses including the 
radiative corrections from the effective Higgs potential are given in table 4. 
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• k,(M^)=0.1 




Figure 5: The lightest Higgs mass including one- loop radiative corrections from quark-squark 
and lepton-slepton loops. Solid (dash) lines are computed by taking various values of ki (see 
legend), Mq = 3 TeV and Ma = 90 GeV. 



3.4 Estimation of the lightest Higgs scalar mass from effective 
potential 



It is well known that radiative corrections coming from the one-loop effective potential can 
give an important contribution to the Higgs masses and push the lightest Higgs above the 



mz threshold |jT8[. These corrections for the case of the MSSM plus a singlet (NMSSM) have 



also been calculated previously||T9[]. In this section we give a quahtative estimation of the 



lightest Higgs mass. This calculation will be very similar to similar calculations in NMSSM. 
The dominant corrections come from the top and stop loops. We calculate a bound in the 
case of large ii approximation. We use the fact that the lightest Higgs mass must be less or 



equal to the minimum of one-loop improved CP-even Higgs mass matrix in the basis (5*1, 6*2). 
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In other words 



9 9 / 9 kfv"^ .9„^x 3, 9 9.9^/, m'}^'m% ( At + li cot m"}, 
m|i < m| cos^ 2f3+^ sin^ 2/5 +7r^/i>t sin^ In + \* ^ . , / In ^ + 



(21) 

Consequently for the spectrum given in table || we get 

m5i<120GeV. (22) 

We have checked that upper bound is saturated when the mixing between the doublets and 
the singlet is negligible that is for small values of the coupling ki. 

In table ^ the numerical values of the one- loop corrected Higgs masses including top-stops, 
bottom-sbottoms and tau-staus effects are given when ki{Mx) = and ki{Mx) = 0.1. 
Lowering the value of the Yukawa coupling ki leaves the rest of the spectrum (other than the 
Higgs scalars) practically untouched. The values given in table ^ are for the case Mq = 0. 
See figure 5 for the case of Mq = 1 TeV. Very large values of Mq will imply a general increase 
in the supersymmetry spectrum beyond the TeV range and may not be experimentally 
interesting. This in principle can also slightly raise the value of msi via the logarithmic 



dependence on stop masses in (pT]) which can be easily estimated. We skip it in this paper. 
However on the brighter side this is not the case if we take ki to be large. This can be seen 
from figure ^ where the values with ki(Mx) = are also plotted. In the large ki case the 
lightest Higgs mass decreases with Mq. As we increase Mq the values of tan/5 < 7 becomes 
incompatible with the current experimental lower bound on m^i. 
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Figure 6: Taking M2 = 90 GeV, we have plotted the parameter space which gives rise to 
charge and color conserving minima (unshaded region) in the Mq — tan (3 plane. When all 
the Yukawa couplings are at the fixed point tan/5 is typically mt/mb ~ 40 — 60. Note that 
there is no charge and color conserving solutions in this case. In this figure we have taken 
kl{Mx) = 47r. 



3.5 EWSB without Charge and Color breaking 

Stronger constraint on the parameter space at the unification scale are obtained demand- 
ing that no charge and color breaking (CCB) minima of the scalar potential is present at 



low energy [0]. Note that we cannot have charge and color breaking minima whereas we 
would like to have radiative EWSB and thereby a correct mz- Large Yukawa couplings 
drive the soft supersymmetry breaking mass terms to negative values provided the super- 
particles whose masses are being considered interacts via the large Yukawa coupling under 
consideration. For example a large top Yukawa coupling drives the Higgs masses to negative 
values triggering radiative electroweak breaking. From(^) we see that ESSM has a number 
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Figure 7: Taking M2 = 90 GeV, we have plotted the parameter space which gives rise to 
charge and color conserving minima (unshaded region) in the Mq — tan/5 plane. Note that 
most of the plane is excluded. Region (b) is due to the effect of a vector-like slepton where 
as (a) is due to three light generations only. In this figure we have taken ki{Mx) = 0.1. 



of Yukawa couplings at the unification scale which couples to the colored superpartners. 
Hence there is a possibility that their large renormalization effects may lead to masses of 
superparticles with weak isospin and color quantum numbers which lead to charge and color 
breaking. This will be the situation when Mq and/or tan/3 are large. Note that large tan/? 
region has a large bottom quark coupling. To get the constraint in the plane Mo-tan /5 we 
verify that for each set of parameters corresponding to a correct EWSB minimum there is 
no deeper charge or color breaking minimum. The allowed parameter space obtained when 
ki{Mx) = is given in figure We have also included the constraint rrisi > 95 GeV 
which forbids low values of tan j3 whereas the CCB constraint excludes large values of tan (3. 
The constraint on the low values of tan/? can be evaded if we reduce /ci(Mx) (figure |^) . 
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Reducing ki implies allowing larger values of the couplings Zf and z'^. This turns out to 
be enough to affect the masses of the extra generation sleptons which are the eigenvalues of 
the mass matrix 



time being in the context of electric charge breaking. Color breaking can be understood 
very similarly. For a vector- like scale of Mi ~ Mq/3 ~ 1 TeV the off-diagonal terms in 
the above mass matrix can compete with the diagonal ones. This increases the splitting 
between the eigenvalues and even rendering one of the eigenvalues negative. Coming back 
to RGE, this would be the situation for most of the parameter space when ki{Mx) = 0.1 
as shown in figure |^. Squared mass parameter of a neutral slepton becoming negative may 
only indicate that it would also get a VEV, and it may affect the values of the other VEVs. 
But a negative charged slepton squared mass parameter would imply a new source of electric 
charge breaking minima in ESSM. This should be understood parallel to EWSB. Note that 
a negative squared mass parameter of the doublet Higgs mav not be enough to generate a 
EWSB minimum. 

Thus in the small ki case save the parameter space lost due to CCB a small region 
corresponding to large values of tan (3 is all that remains. We point out that there is no such 
constraint due to the vector-like spectrum for larger values of ki. 




where m|;. are the running mass parameters. Let us discuss the slepton masses for the 



26 



4 Conclusion 



We question the commonly accepted notion of a unified gauge coupling ax ~ 0.04. If the 
Minimal Supersymmetric Standard Model (MSSM) is extended by including two vector-like 
families (ESSM) the couplings grow stronger than the low energy ones due to the renormal- 
ization effects of the extra matter and unify at a semi-perturbative scale of around 0.2. This 
is actually the only extension of MSSM containing complete families of quarks and leptons 
that is permitted by measurements of the oblique electroweak parameters on one hand and 
renormalization group analysis on the other. The former restricts one to add only vector-like 
families whereas the latter states that no more than one pair of families can be added to 
maintain the perturbative unitarity up to the unification scale. In ESSM the weak SU(2) 
couphng grows by a factor of six at the unification scale compared to the weak-scale value (fig- 
ure |l].a). The four dimensional string coupling may have a similar intermediate value which 
is large enough to make the dilaton stable as was conjectured by the previous authors[l|]. 
ESSM has a unique pattern of the Yukawa matrices which is motivated by preon theories. 
The vector-like matter and normal matter has off diagonal Yukawa couplings whereas the 
normal three generations do not have Yukawa couplings among them at all. This leads to a 
see-saw like picture of the fermion masses. Below the mass scale of the vector-like generation 
a hierarchical mass pattern of chiral fermions emerge. If we fix all the Yukawa couplings to 
be large at the unification scale we get unique predictions of the low energy fermion masses 
when the Yukawa couplings approach their "quasi-infrared fixed points" at the scale of the 
top quark mass. On the contrary the renormalization effects of these relatively large Yukawa 
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couplings have non-trivial effects on the unification of gauge couplings. Keeping this in mind 
we have also performed the renormalization group evolution of the gauge couplings taking 
into account the Yukawa effects at the two-loops. If we assume the universality of the soft 
supersymmetry-breaking parameters at the unification scale renormalization group evolu- 
tion enable us to determine the supersymmetry spectrum at low energy quite easily. Note 
that due to the presence of the heavy generations the renormalization of the superparticle 
mass parameters are considerably different from that of MSSM as we would expect. This 
makes ESSM distinct from MSSM from the point of view of collider searches. The first and 
second generation squarks do not have large Yukawa renormalization hence they experience 
pronounced QCD renormalization which make them heavy (figure 3.b). 

A further question will be to get the correct radiative electroweak breaking. We point out 
that the mass of the vector-like generations is actually linked to the electroweak symmetry 
breaking mechanism by the approximate relation Mq ~ (zq/ki)/!. An electroweak symmetry 
breaking minimum which fits the mass of the Z boson exactly cannot be obtained in the 
"quasi-infrared fixed-point" scenario of the Yukawa couplings if we like Mq to be below 100 
TeVs. This is too large to be interesting experimentally. Thus at least some of the Yukawa 
couplings must have smaller values. The fixed-point scenario naturally have a large tan/3. 
We show that large tan/3 region suffers from the presence of Charge and Color breaking 
minima for any value of the vector-like scale from mz up to Mx- So also in this case we find 
that to get a global charge and color conserving minima we must give up the assumption of 
all Yukawa couphngs at their "quasi-infrared fixed point" in the case of ESSM. 
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5 Appendix 

5.1 RGE coefficients of Yukawa couplings 

The RGE for the Yukawa couphngs (at any order) are given by the following general expres- 
sion which is perturbatively exact 

16vr'^ = htilL)] + {ir)\K, + K,1H^ ■ (24) 
In (1^) is the Yukawa coupling for the right-handed field the left-handed field Lj and 



the scalar Hm- The functions 7^ are the anomalous dimensions for the superfield A||20|. We 
define /i^ as the Yukawa coupling matrix when all the superfields enter in the vertex and 
hjj when they leave. 



We have used the two-loop anomalous dimensions to evaluate the expression ( P^ ) which 
can be splitted as follows 

(74 = (7i)} + Y^(7i)5. (25) 
The component one-loop and two-loop parts are given by the expressions 

= h^ht; -2glC2{Ra)5.,, (26) 
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ill)) = - {Kkil^B^fiK' + Kk^lKi) + 1haC,{Ra)g% - 2{^^l^))C,{R^)gl . (27) 

We have denoted Qa as the gauge coupUng ha as in the one- loop gauge /9-function and C2{Ra) 
as the quadratic Casimir operator for the Ra dimensional irreducible representation. The 
anomalous dimensions can be expanded in terms of four Yukawa coupling matrices related 
to equal number of Higgs bosons (at the unification scale they can be thought of a 10 and 
two singlets of SO (10)). The index m assumes the values m — Hi, H2, Hs, Hx whereas 
f = u, d, I, u. The up sector Yukawa matrix can be re-expanded in terms of the individual 
matrices which are 

/ \ 

c = , (28) 

V / 

/ 

I , (29) 

V y'/V2 





/if^ = I Zu \ . (30) 
z'jV2 

We will have similar expressions for the down-quark sector and similarly for the leptonic 
sector. The normalization factor y/2 avoids over counting when we sum over the index 
f ^ u, d or f ^ I, u. 

The one-loop anomalous dimensions for the quark sector are as follows 



i^lh = 2(/.f/.f^),, + (/.f/^f^),, + (/.f/.f^),,-(^^,? + ^^32)5,,, (32) 
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+{h^''h^% - i^gl + \9l + l9^)8^, , (33) 

7i = (/^^'C^)(3,3) + (C^'/^f)(3,3) - - (35) 

7i = (/^^/^^')(3,3) + (/^^/^^')(3,3) + {h^'h^^'\,,^) + (/..^^/.f ^3,3) 



^ ( <t' 4 ) ' 

The two-loop contributions to the anomalous dimensions for the quark sector are as follows 
(7?).. = -(2/^f(7^„+7k)e'+^f^(7^„ + 7i.J/^f^V/i^(7^„+7K)^^')., 



- [-^91 + 3^?3 j {ll\3 + (^Y^^i^r + 3&3^3 j ^ij 

(7|)., = -(2/if(7^, + 7^J/^f'+/^f (7^, + 7L)^f'+^^(7^, + 7k)^^0, 
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.30"' 2"" V30 2 3 



(3,3) 



+ 3^1 j 7^ + Y^Mt + 
7!) = -(/^^'(7U7k)^?^ + C^'(7i5 + 7l.J/^?^ 



(3,3) 



15^^ + 3^3 j lo + -,^19. + 3^3 



7j = - (^^(7^„ + 7k)^^' + ^^(7^. + 7k)^?^ + h^'i^L + ^hsK'' + h^'hL + 7^ 

" IsO^^ + 2'' + 3^^j ^« + 30^^^^ + 2^''' + 3^^^^- 

Two-loop anomalous dimensions of leptons can be obtained from the above expressions with 
the replacements u ^ u, d ^ e and q ^ I. 

We have set the Yukawa couplings for the first and second generation to zero. They can 
be included by the replacement of the numbers x/, x'^, yj and y'^ by corresponding three 
dimensional vectors. In this case 3x3 Yukawa matrices given in equations (|28| - |30D will 
become 5x5 matrices. 

The one-loop and two-loop anomalous dimensions for the Higgs scalars are as follows 
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(a) one- loop anomalous dimensions 

(a) two-loop contributions to the anomalous dimensions 

7k = -3Tr[/.f 7i>^^^ + hf^hh^^] - Tr[h^^^lh^^' + /if ^^^/if ^ 

-^i(7k + 7k)^i - {^9l + ^^2) 7k + ^b^at + ^^252' , 
7k = -3Tr[e7i„e' + eV^C] - Tr[h^'llh^^' + ^^7^^] 

-^i(7k + 7k)^i - (^^1 + ^ai) 7k + ^^i^f + ^^2^2' , 

7k = -6Tr[/.f-7i>f + h^/llh^a' + ^^7^^^' + ^%^^^] 
-hlRgki, 

-^37k^3- 
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5.2 Two- loop formulae of threshold corrections in step- function 
approximation. 

(k) 

The two-loop coefficient for the RGE of the gauge couphng bij and including the threshold 
corrections are listed here. 9i is the step function. The mass thresholds are denoted by the 
index i = w,g,qL, ur, dR, II, cr, h, H and similarly the vector-like thresholds i — Q, L. 
Here Uf is the number of chiral generations and ny is the number of extra generations, 
(a) Gauge contribution 



J21 



■>22 



J31 



-'32 



'19 



9 



1 



■Onr 



+ 



9 

50 



32 ^ 
75 



bi2 = 



(81 137 
T00^^ + 300^«, 

1 + 9h+0j^{1-9, 



+ 



9 



20 ^^"^ 



10 

+nv (3 - 9^) (^^9l + ^9q 

/AA / A 



.)) 



(— + (2-9.)( —9., + —93 + —9u 
V15 ^ ^' V15 15 15 
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+ny (3 - 9g) -9q 



— V 15 ^' 

+4(l + 9H + <'i.(l-9»)) 



+».(3-fc)(|j0. + i0«) 



&23 = 



136 64 
/13 

+ny (25- lift, 

^/ (4 + 4(2 
+nv4 (3 



4 /49 /I 1 

-0^ + ^/ (y + (26 - 33^^) [^-9i^ + 



6 2 

1 , 3 



(30 + <2 



30 

% + ^ (2 - ^.~.) 



/I 1 4 > 

'^^ ( 30^^"- + 15^'^« + I5S 



(54) 



(55) 



(56) 



(57) 



(58) 
(59) 

(60) 
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+nv (40 - Oq. 



(61) 



(62) 



(b) Now the coefficients : Let us become careful here. We define 7j to be the anomalous 
dimensions without the gauge contributions which are as follows 



E«2 

k 

E«3 



6 /4 



5 V3 



+ 2 (^k +7^ 
3 (Tr7i + 7^) + Trf, + 7i + + 7),, , 

Tr7i + Tr7i- + 7I, + 7I, + 2 {Tr% + 7J) . 



(63) 
(64) 

(65) 



The vector-like superfields are massive at the scales (Mg, Mj)j. After the rotation of the 
Yukawa matrices we get 







( htifi) 



















)■ 


(66) 


















= /if(/i<MQ; 


= 


0, 


(67) 






( hifx) 









/if (/x<Mq) 










)■ 


(68) 
















/if (/X<Mq) 




= 


0, 


(69) 






( Kifx) 









/if (/i<M^) 










)■ 


(70) 






[00 









/if (/x<Mi) 


= hf-{fI<ML) 




0, 


(71) 



with our boundary conditions (|^). 
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niiiGeV) 


tan /9 = 5 


tan /9 = 30 


tan /? = 50 


Mo, Ml/2 


0, 1115 


0, 1083 


0, 952 


Ao 


-2768 


-1922 


-1800 


V\ 


4828 


4964 


6389 


/i = ki{mz)vx 


1169 


1091 


953 


B = k2{mz)vx + Ai{mz) 


297 


34 


4 




64 


73 


81 




1342 , 1425 


1057, 1490 


439, 1171 




1351, 958 


1051, 1458 


436, 1455 




1345 


1057 


441 


rfiq., i = 1,2 


1185 


1157 


1042 


m^., z = 1,2 


536 


522 


466 


rhr^, i = 1,2 


355 


345 


305 


rht^, rht^ 


653, 974 


721, 879 


680, 750 




913, 1074 


840, 984 


617, 816 




634 


610 


538 


rhr-,, rhr^ 


355, 645 


265, 624 


204, 568 


Xl A.2 


87, 1175 


89, 1097 


89, 960 


56, 90 


58, 90 


58, 92 




957, 957 


1094, 1094 


1168, 1174 




1530 


1710 


1441 


M3 


357 


355 


348 



Table 3: Supersymmetry mass spectrum in the ESSM for our case. Mq = 3 TeV and 
M2 = 90 GeV. The values of the mass parameters at the unification scale and the values 
at the 100 GeV scale are included. The Higgs masses are tree-level masses. Correct EWSB 
is guaranteed for this mass spectrum. Corresponding Yukawa couplings are given in table 
2. The lightest Higgs mass m^^ goes above the experimental lower bound of 95 GeV when 
one-loop radiative corrections are included (table 0). 





ki{Mx) = 3.5 


k,{Mx) = 0.1 


mi{GeV) 


tan /3 = 5 


tan P = 30 


tan (3 = 50 


tan (3 = 5 


tan/3 = 30 


tan (3 = 50 




120 
1061 , 1423 
1090, 934 

1061 


105 
1227, 1490 
1219, 1461 

1190 


107 
961, 1171 
959, 1455 

799 


115 
1342, 3950 
1342, 6840 

1344 


119 
1176, 2124 
1176, 3678 

1147 


116 
847, 1439 
847, 2493 

698 



Table 4: Higgs masses after including one-loop radiative corrections to the scalar potential. 
Mq = 3 TeV (Mo = 0), M2 = 90 GeV (experimental lower bound on the wino mass). 
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